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^^ ■ Abstract 

Cn| . 

This paper is devoted to the long-time asymptotic analysis of the Cauchy problem for the modified 

rS ,■ Korteweg - de Vries equation with a step-like initial function, which rapidly tends to different constants 

as X — >■ ±cxj. 

First investigations of this problem were done by E. Khruslov and V. Kotlyarov in 1989. By using 
the technique, developed by E. Khruslov for the Korteweg - de Vries equation in 1976, they studied 
the long-time asymptotic behavior of this problem solution in the domain x > 4c^t — const log t, and 
T^ I found (discovered) that the solution breaks (splits) up into a train of so-called asymptotic solitons in 

♦rp^. the subdomain 4c^t — const log t < x < 4c^t. 

r^ I Further, in 2010, V. Kotlyarov and A. Minakov studied the behavior of this problem solution 

"ti . in the whole x, t plane with the exception of small neighborhoods of the leading and trailing edges 

\x + 6c^t\ < et and —et < x — Ac^t < 0. It was found that the solution of this problem in the domain 
(— 6c^ -I- e)t < X < [Ac^ — e)t is described by the modulated elliptic wave, as was earlier predicted by 
A. Gurevich and L. Pitaevskii in 1973. 

In this paper the long-time asymptotic behavior of this modulated elliptic wave is studied in the 
JZ^ ' small neighborhood of the leading edge 4c^t — const log t < x < ic^t. We show, that this modulated 

^— V . elliptic wave also breaks up into the train of asymptotic solitons, which are similar to Khruslov solitons, 

J — ' but differ from them in phase. 

Key words: nonlinear evolution equations, modulated elliptic wave, asymptotic solitons. 
•^ ' Mathematical Subject Classification 2000: 35Q15, 35B40. 
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1 Introduction 

The history of long-time asymptotic analysis of the Cauchy problem for the modified Korteweg - de Vries 
equation 
^.' qt{x,t)+6q^{x,t)qa:ix,t)+qxxx{x,t) =0 (1.1) 

with initial function of step-like type, i.e. 

I as x — )• -|-oo, 
q{x,0) = qo{x) ^ { (1.2) 

I c as X — )• — oo 

is quite long and goes along with long-time asymptotic analysis of the Cauchy problem for the more 
famous Korteweg - de Vries equation. Physicists have understood qualitative description of the solution 
since pioneer work of A. Gurevich and L. Pitaevsky [1] (1973). Namely, x,t - plane is divided into 
three domains. In the left and the right domains the solution tends to constants, and in the middle 
domain it tends to a modulated elliptic wave. Between these domains there are transition regions. In the 



transition region near the leading edge the chain of asymptotic sohtons run. In this direction R. Bikbaev, 
R. Sharipov, V. Novokshenov actively worked (see [2]-[H]). All these works were done in the framework 
of a heuristic Whitham method. 

There were not any rigorous mathematical papers on this theme except for papers concerning soliton 
domain. It was proved that near the wave front there are arise the so-called asymptotic solitons, which 
are generated by the edge of continuous spectrum, as opposed to usual solitons generated by discrete 
spectrum. For KdV equation it was done by E. Khruslov |15) . |16| (1975, 1976), and for MKdV it was 
done by E. Khruslov, V. Kotlyarov |18) (1989). The review of results in this direction for many other 
integrable equations can be read in |19) . |20) and in the bibliography of this papers. 

From the other hand, for more than 20 years the method of Riemann - Hilbert problem and the 
steepest descent method ([?]) have been being actively developed. Recently these methods were applied 
to study solutions in other domains of x,t half-plane, not only in soliton domain ([l3], |9] - |12) . |21j - 

Let us consider, in particular, the Cauchy problem for the modified Korteweg - de Vries equation 

qtix,t) +6q'^{x,t)q^{x,t) +qxxx{x,t) = 
with initial data, which rapidly tends to and some constant c > as x — )• ±00 : 

I as x — > +00, 
qo{x) = q{x,0) -^ < 

\c as X — ;■ —00. 

Suppose that there exists a solution of this problem q{x, t), which rapidly tends to initial constants and 
c as X — 7- ±00 for all i > 0. We interest in long-time asymptotic behavior of the solution of this problem. 

This problem was considered earlier by E. Khruslov and V. Kotlyarov in 1986 |18) . and also by 
R. Bikbaev [6] in 1995. 

In |18) the authors extended Khruslov method, which was originally developed for the Korteweg - de 
Vries equation, to the modified Korteweg - de Vries equation. They studied the asymptotic behavior of 
this problem solution in the neighborhood x > 4c^(t — const log t) and calculated explicit formulas for 
asymptotic solitons, which run in the domain 4c^(t — constlogt)x < AcPt. Just this result is compared 
with our result in this paper. 

Also this problem, with more general type of initial data, was considered by R. Bikbaev in a short 
note [6] in 1995. He used a heuristic Whitham method and his results are based on an analysis of the 
Whitham equations and the theory of analytic functions on a hyperelliptic surface. 

Later this problem was revisited by V. Kotlyarov and A. Minakov |17) in 2010. There the nonlinear 
steepest descent method for oscillatory Riemann - Hilbert problem was used and applied for rigorous 
establishment of the asymptotics. In this paper the asymptotics were studied in the whole x, t plane with 
the exception of the neighborhoods of the leading and trailing edges, that is, in the domain x < — (6c^+e)t, 
-{6c^ -e)t<x< {Ac^ - £)t, X > AcH. 

In particular, in the domain — (6c^ — e)t < x < (4c^ — e)i the asymptotics of the initial Cauchy problem 
solution is described by a modulated elliptic wave: 

X / C £ C £ \ 

Theorem (Kotlyarov, Minakov) Let ^ = . Then for ^G I 1 , J the solution of the 

IB V problem M.l\) . M.2\) takes form of a modulated elliptic wave 

q{x,t) =qei{x,t)+o{l) , t -^ 00, 

where 

-e (^i + itBgifiO) + iA(/(0)|r(/(e))) , X 



qei[x,t)-Vc / (?) e(iti?,(/(0) + iA(/(0)|r(/(6)) ' ^ ~ I2t ^'''^ 

Here Bg(d),T{d), A(d) are defined by \2.l4^ )- [2.14\ ) and the function /(?) is defined by lemma \KT[ It 
worth to mention that the point \d = i/(?) is an analogue of a branch point of an Whitham zone, r - is 



the b - period of the corresponding Riemannian surface, Bg - is the b - period of a Abehan integral of 
the second kind, and A - is an integral on a Riemannian surface. 

Here e is arbitrarily small positive number. It was introduced for a technical reasons. The question is 
if it just a technical restriction or more significant limitation, and there must be another formula in the 
domain (4c^ — e)t < x < x. 

Notice that the explicit formula, which determines qei{x, t), is valid up to the leading and trailing edge, 
i.e. in the domain —6c^t < x < A(?t. This suggests to study the behavior of q(.i{x,t) in the neighborhood 
of the leading edge x = Ac^t and to compare the result with one obtained in [?]. 

In this paper the long-time asymptotic behavior of this modulated elliptic wave ge;(x,t) is studied in 
the small neighborhood of the leading edge Ac^t — constlogt < x < Ac^t. The main result of this paper 
is: 

Theorem 1.1. Let N >1 he any integer number, e > is a fixed small number. Then as t —)• oo 
A uniformly for t and x such that 

AcH log t < X < Ach - log t 

2c e. - - 2c ^ 

the following is true: 

qei{x,t) = g ^^^^ ^2c(x - AcH) + (2n - i) logt - a„(rE,t)) ^ ^ UJ ' ^^'^^ 



B uniformly for t and x such that 



Ac^t logt < X < Ac^t 

2c 



the following is true: 

qei{x,t) = 0[t 



Here an{x,t) is an auxiliary well-defined function, which has the following asymptotic behavior as t —t- oo : 



where 

v= 1 



Ac'^t' 

For any 7, 5 such that < S < j, 
an{x,t) is bounded in the interval A(?t — 7logi < x < AcP't — 5\ogt. 

In the section m we discuss this formula and compare it with the formula for q[x, t) from |18j . 



2 Preliminaries 

Let us consider the initial-value problem 

qt{x,t) +6q'^{x,t)qx{x,t) +qxxx{x,t) = 0, 

I 0, X — )• +00, 
q[x,0) =qo{x) -^ < 

c, X — ;• —00, c> u, 



with initial function which rapidly tends to their limits, namely 

+00 

/ \qo{x) - C(|e2|^l('+^)d2; + / \qo{x) - cje^^^'+'^Mx < 00. (2.6) 

-00 

Here / satisfies the inequality I > c > 0, and cr > is an arbitrary positive number. We suppose that 
there exists the solution of this initial value-problem and that this solution converges to its limits with 
the first moments: for any t 

+00 

/N|,(.,0-c|d.+ /.k(.,t)|d.<oo. (2.7) 

-00 

Let a~^{k) be the standard transmission coefficient (|17)). It has the following properties: 

1. a{k) is analytic in the upper half-plane C+ with the cut [0,ic], extended continuously up to the 
boundary with the exception of the point ic, where a{k) may has a root singularity of the order 
(fcTic)-^/^ 

2. a{k) satisfies the symmetry conditions 



a{—k) = a{k); 

3. V/cG (ic,0] a±{k) / 0. 
We suppose the absence of solitons, i.e. that 

VA:GC+\[ic,0] a(A;) / 
and that a{k) has a root singularity at the point ic of exactly order (fc — \c)~^''^: 



The asymptotic behavior of this Cauchy problem (jl.ip . ()1.2p solution was studied in |17| . The x, t- half- 
plane was divided into 3 regions with qualitatively different asymptotic behavior of the solution, namely 
X < (— 6c^ — e)t, (— 6c^ -|- e)t < x < {Ac? — e)t, x > Ac^t, where e > is arbitrarily small positive number. 

The most interesting in asymptotic analysis is the domain (— 6c^ +e)t < x < (4c^ — e)t. In this region 
the Cauchy problem solution is described by a modulated elliptic wave, according to formula (|1.3p . 

Our notations 
• 0-function Q{z\t). 

In the sequel we use 0-function in the next form: 
for any r < and z £ C 



Q{z\t) = y^ exp < -rm + zm > 

TTl — — rv~i ^ ' 



(2.8) 



m=— 00 



0-function has the next property: 
for any integer n, I 

e(z + 27rin + tI\t) = Q{z\t) exp i --rf - zl\ . (2.9) 



• Real parameters d and fi. 

Function d = f{^) : — ^, % 



[0, c] is defined from the following system: 
1 ; 



(^^-AV)J^i-^dA = 0, 



(2.10) 



Lemma 2.1. For $^ G 



-c2 c2^ 



C' , 2 '^' 

i/iere exists one to one, increasing and continuous map f : 



(2.11) 



-c2 c2 



[0, c] which defines the unique solution of equations i2.10\) and h2.11\) by formulas: d = /(^) and ^ 

is also continuous and increasing, and their boundary values are 



'V3 



0, /i(0) = 0, m(c) = -^. 



^(/(O)- The map fi : [0, c] H 

as follows: f [^j = ^' f [~ 

The proof of lemma [2?T] is given in Appendixes. 

• Functions Bg{d),T{d), A{d). 

For d G (0, c) let w{k, d) = y^(/c^ + c^){kP' + d^) be an analytic in k function in the domain 

k E C\ ([ic, id] U [— ic, —id]) witch is fixed by the condition w(0, d) = cd. Then we define Bg{d) , T{d) , A{d) 

as follows: 



Bg{d)=2AJ 
id 



{k'^ + i?{d)){e + d'^)dk 

w+(A;,d) 

id 



r{d) 



-vri 



dA; 



w+(fc, d) \ J w{k,d) 
id \0 



dA; 



(2.12) 



(2.13) 



A(d) 



id 



log {a^{k)a-{k)) dk 
w_|_(A;, d) 



dk 



w{k, d) 



(2.14) 



3 Proof of the Theorem 11.11 



Theta- functions in ()1.3p are slowly-convergent as ^ tends to — because in this case t(/(^)) tends to 

(see lemma ED in appendix). To overcome these difficulties, i.e. to make theta series well converged, we 
can use Poisson summation formula and rewrite ()1.3p in term of rapidly convergent theta series. Indeed, 
we have 



G(z|t) = G 



f27riz 



V ''" 



47r^ 




2tt I -z 



Then the model approximation (|1.3p of the solution takes the form: 
qei{mi,t) = 



(3.15) 
(3.16) 



V^^^7^exp(ziMi))^ilM)(.(,e) + r^"v 2 



0,I%P)(^(i,^) + l) r*{f{0) 



Qi^^HP^z{t,Or*{f{0) 



1 47r2 

where z{t,^) = - {tB„{f{^)) + A(/(0)) and T*{d) = — — . We emphasize that now r*(/(0) ^ -oo as 
vr ''"(") 

^ — 7- — and as a result theta series become well-converged, 
o 
Let us introduce a function F(.,.). For any r* < and z G C we set 

F(r*,z) = Vc^-/^2(r*)exp^ + -(^ + l) ^ ^. (3.17) 

Here h{.) is defined by lemma \5A\ Thus we can rewrite the formula (|3.16p by using the function F{., .) : 

q,l{x,t) = F{T*{f{0),4t,0)- 
Lemma 3.1. For any integer N > 1 the following formulas is true uniformly in z as t* ^ — oo.' 

Fir*,z) = Y: ^ r rH-i-^+2n) + Q {^^) ^ 0<z<2N, 

n=l cosh ^ i ^ ^ 

F(T*,^) = 0(e*), =^<^<\- 



Proof. We have that 



e I -•(' + ')|r- 



e I ^1.- 



Taking into account the definition of function (see formula ()2.14p ). we can easily verify that for r* — t- — oo 

1 5 

the next formula is true uniformly for < z < — : 

■^ 2 - - 2 

,, , / / .-xx .♦(.+!) 0(e^*(i-^))+e^+0(l) 
F(t\z) = 4c ( 1 + O ( e- ) ) e^^ ^ Jj--, ) ,,,^ , 

^ ^ ^^ 0(e^*(2-))+e^^l^ + l + 0(e^^i^ 

Simple calculations lead us to the next formula: 

4c (l + O (e^)) (l + O (e^ + e^^^ 

-^ V 5 ^j ~ r*(l-z) -r*(l-z) / T*(7-3^ T*(l + 3z) 

e i +e 4 +0e 4 +e * 



and then 

4c (l + O ('eT^')') A + O (e^ + e^-^ 

^V ) ^) = -- —7 / T*(7-3z) T*(l+3z) 

Now let us note that 

T*(7-3z) T*(l+3z) 

^ , e 4 +e 4 \ / ^"(3-^) r*(l+z) 

01 —^^ l=0(e 2 +e 2 

cosh — 5^1 — - 



Then 

F(t,2;) = -7- — ^ 1 + e4+e2 +e 2 +e 4 +e 2 

and 

. * s 2c / ( T* r*z T*(2-z) 

^(^ ' ^) = ?^?T3IT 1 + O e- + e— + e^^ 

cosh — ^-^ — '- \ \ 

Therefore 

(r^ 2^ r*(2-z) 

e 4 +e 2 +e 2 
cosh ^ ^ 

The argument of the 0-estimate in the last formula may be change with: 

e 4 tor < z < — , 

2 - - 2' 

T*(2-2) 1 

e 4 , for - < 2; < 1, 

2 - - ' 

e^ , for 1 < z < -, 

- - 2' 

^*(3-') ^ 3 ^ ^5 

e 4 tor - < z < -. 

2 - - 2 



Therefore 



and 



F(r*, z) = ^^ + O (e^) for < z < 2, 

cosh "^ \^ ^ ^ ^ 

F(r*, z) = O fe'r') for ^^ < z < -. 

2c / T* \ 

Also we note that — — ^ = O ( e^^" ) for z G M\(0, 2). 

Besides this, due to property (j2.9p of G-function, we obtain that F{t* , z) = F{t*, z — 2n) for any integer 
n. Therefore 

F(t*,z) = -^^ TTTT + Ofe'r') for 2n - 2 < z < 2n, 

coshi^tiz£±M V ; - - > 

and 

2c / * \ 

O(eT-) for z G M\(2n-2,2n). 



cosh -'(-^7+^") 



So, we obtain that for any integer A^ > 1 



F(r-,z) = ^^ - ^ _^ +0(e^), < z < 2Ar. 

n=l cosh ^ ^ ^ ^ 

F(r*,z) = 0(e^), =^ < ^ < \- 



n 



Now we have to express in terms of variables x, t the intervals z G [0, 2A^] and z G [—2' 2]' To do this 
let us define small parameters r/ and v. 

rj = r](d) = 1 , and v = v(S,) = 1 ^, 

c c^ 

which tend to zero as .^ — )■ c^ /3. The inverses of these functions are 

d{v) = cil-7j), and C(^;) = ^!il_^. 

In the appendixes we find the connection between v and r/ (see the formulas (|5.4U|) . (|5.42p ). and make 
asymptotic expansions of Bg{d{r])), A{d{r])), T*[d{r])) in rj (see the formulas ()2.12p . ()5.34p . ()5.30p ). 

Lemma 3.2. 

, ^, ,, ^cHv 1 ^ / 1 tuloglogA 

logi 2 \^log- log V J 

V 

Proof. Recall that 

/ . i/r. ^^^ tB Jdir])) + Mdir])) . 1 21og2 ^/ o, 1 1 

z{t, r\d{7]))) = ^^ ^'^' ^^^ = 8chr] - - + V + O ( tif log - + 



VT ^ log- ^ '^ ^°§ "? 

?7 

From this formula we can deduce that for any positive t and sufficiently small r] : z > . Further we 

will use this remark. 

By substituting in the last formula for z expression for 77 in u we obtain the equation (j3.18p . D 

Correspondence of intervals in z and in x,t. 

Lemma 3.3. 

A. For any integer N > 1 and any sufficiently small positive e there exists T = T{e,N) such that for 
any x, t which satisfy the conditions t > T and 

Ac^t ^ \o^t<x<4:C^t ^logt 3.19 

2c 2c 

the following formula is true: 

0< z{t,i) < 2N. 

B. For any integer N > 1 and any sufficiently small positive e there exists T = T{e, N) such that for 
any x, t which satisfy the conditions t > T and 

Ac^t — \ogt<x< Ac^t (3.20) 

2c 

the following formula is true: 

-\<z{t,0<'^N + \- 
Proof. Let 7 and 5 be a some positive numbers such that 5 < 7 and let us consider interval 

4c^t-7logt < 2; < 4c^t-(51ogt (3.21) 

or 

7logt X ^ (51ogt 
1 ■ — 2_ < < 1 2_ 

4c2i 4c2i 4c2i 



X 

We remember that v = 1 tt- and then the last formula can be rewrite in the next form: 

4c^t 

We see that for x, t, which lay in the interval 

Ac^t - 7 log t < X < Ac^t - 6 log t, 

if t — 7- oo then t; — )■ 0. Further we will use this fact. 

For t and v such that 

5logt 7logt 

11 1 

< T- < 



also 



and then 



and 



logt loglogt-log^ ,1 logt loglogt-log^ ' 

^ logt ^'^ V ^ BP 

2cS 8chv 2c7 

< T r < 



r log t-log ^ loe - . log log t-log %2. 



log* 



Lf# logi; 1 



2c5 8c^tv 2c7 8c^tu 

+ z < z < -^ — i- + z. 



_ log log t-log :^ log i log log t-log ^ log i 

logt -•- logt 

And then, due to formula ()3.18p . we obtain that for any positive e there exists T = T{e,'y,6) such that 

for any t > T, for any v which satisfy p.22p . and uniformly for 5, 7 from the interval — < 5, 7 < 

Ac 2c 

the following inequality holds true: 

2c6 e < z <2cy he 

2 - - ' 2 

A.T 1 ,. ^ + 5 2iV + i-e 

A. Wow let us take = —^ — , 7 



2c ^^^ 2c ^^ 

If we take into account that (|3.22p is equivalent to (|3.2ip . then we obtain that for any positive e there 

exists T = T{e, N) so that for any t > T and any x such that 

, 2iV + i-e , i + e 

4c^t log t < a; < AcH - log t 

2c B - - 2c ^ 

the following is true: 

< z < 2iV. 

B. If we take now 7 = and take into account that always for big t : z > ^, then we obtain that 

for any positive e there exists T = T(e) so that for any t > T and any x such that 

Ac^t logt < X < Ac^t 

2c 

the following is true: 

-1 1 

— < z< -. 

2 - - 2 

D 



Evaluating of cosh argument in the formulas for q^i 

Define a„(x,t) through the relation 

'^ ^^~ " =2c(x-4c2i)+ ( 2n- - ) logt-a„(x,t) 



4 ' ' V 2 

Then we have the following lemma 
Lemma 3.4. 

an{x,t) is bounded for Ac^t — 7logt < x < 4c^i — 5logt, where < 6 < ^ is arbitrary numbers. 

T*(z — 2n + 1) 
Proof. Let us express in terms of v,t and therefore x,t. We know that 

8 , , / log log i + 3 log 2 log log i + log 8e /log^logi^ 

r* = -41og - + Or/ = 41ogt; 1 - ^ — + ,% ^- + O f , " 

r? \ log I) log V \ log'^t' 

tBa + A 



z 



vr 



1r q 3 , n^ 2i ^ 8c^^ / log log -+ log Be \ /ulog log-^ 

-Sg = Sc^r/ + O (?7 log??) = ^ 1 + f +0 -^-^ 

vr ^ logi \^ logr; J y log'^u ^ 

Ia = ^ - ^i^ + O f i^sbgi \ 

TT 2 logU I log^V I 

And then for t,v, which satisfy (|3.23p . in forward, but careful calculations we obtain that 
T*(z-2n + l) o / 1\, / 1\ logi 8c^tv 

/ 7\ /log2logi\ 

+ 6n - - log 2 + O ' ^ ' 



2/ \ logti / 

And by remembering that v = 1 — j% is function of x, i, we obtain that 

T*(z-2n + l) , , 1 ^ ( 1\ , , , 

^ ' - 2c{x - Ach) + 2n - - log t - an{x, t), 



4 ' ' V 2 

where 



IV logi 8chv (^ 7^, ^.^/log^togi 

logf 



„„(,, t) = - (2n - i) log !^ - ^ - (6„ - Z) log 2 + O 



While proving lemma [HTBl we also have proved that an{x^i) is bounded for 4c^t — 7logi < x < 4c^i — Jlogt, 
where < 5 < 7 is arbitrary numbers. D 

By summarizing lemmas l3.H 13.31 and 13.41 we obtain the statements of the Theorem ll.il Q. E. D. 
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4 Comparing of qei{x,t) and q{x,t) in the domain 4c^t — 7logt < a: < 4c^t 

In this section we compare q(,i{x,t) with q{x,t). In |18) the following result was obtained: 

Theorem (Khruslov, Kotlyarov) Let N > 1 and M be such an integer numbers that N 
Then for 

2c 



M + 1 



N 



q{x,t)=2cY,' 



"=i cosh ( 2c(x — 4:C^t) + ( 2n — - ) log i — a„ 



where On = log 



-p(n)-p(n) 
^1 ^3/2 



(n - l)!242n-l(2c)6"-3 p(n-l)p(n-l) 



3/2 



.(fc) 



det[[r(i + j + 6)]L, 



ij=0,k-l ' 



+ t"2 



Here ^o is some constant, which is determined by qo{x). 

We can rephrase the theorem by simplifying expression for a„. Also as both M and A^ are integer we 
can take M = 27V. 

Theorem (Khruslov, Kotlyarov) (rephrased) Let N > 1 be an integer number. Then for 

2, (2iV + l)logt 

X > Ac t 

2c 



N 



q{x,t)=2cY,- 



w=i cosh ( 2c{x — AcH) + ( 2n — - ) log i 



Ctr 



+ 0{t'2 



(4.24) 



where 



din = log 



A2nfitcfn-, nn)T{n+\ 



(4.25) 



1 



„ 2m — 2 Olrn 

Let us consider a curve x = Ac t log t -\ , where m G N. All the summand in the formula 

(j4.24|) become vanishing except for one where n = m. Then for x, t which lie on this curve 



j{x,t) = 2c + 0{t 



-+e 



as i — )• oo. 



Now let us consider q^i{x,t). All the summand in the formula (jl.4p become vanishing except for the 

one, where n = m. According to the formula (|1.5|) . for x = Ac^t log t -\ — and t — >• oo we have 

the following: 



vt 
1^ 



^"""^ /^l + O '^^°s^°si 



8c3 



logt 



11 



vt 1m -\ /^loglogt 



logi ^ 8c3 V log* 

and 

/log^logt 

where 

am = [iva - M log ^^3-^ - Tern - ^ j log 2. (4.26) 

Then 

2c ^ / -i+f\ 2c ^/log^logt 

cosh (am - am (a;, tj) V / cosh (a^- am) V log* 

As Um (|4.26p and 0.^ (|4.25p are not equal, q{x,t) — qei{x,t) does not tend to as t — )• 00. 



Remark. In asymptotical analysis in the domain x < Ac^t (in |17| ) we had to restrict ourselves for 
a smaller domain x < AcPt — et, where e > is arbitrarily small number. It was done for technical 
reasons. We see, however, that the asymptotics from the domain 4c^t — et < x with arbitrarily small e 
and the asymptotics from the domain 4c^t — 7logt < x do not match. This may indicate the presence of 
a transition zone 4c^t — et < x < A(?t — 7 log t. 

5 Appendixes 

4^2 

Lemma 5.1. Let T*{d) = , where r(.) is defined by the formula ()2.13p . Then 

1. ford£ (0,c) :r*(d) < 0; 

2. T*(.) is decreasing function on the interval d G (0,c); 

3. r*(+0) = 0, T*{c - 0) = -00. 

Therefore there exists the inverse ofT*(.), call it h{.) : (— oo,0) — )■ (0, c) : h{T*{d)) = d. 

Proof. Indeed, 

id 

dA; \ I f dk 



-<^)--/^ / 



w+(A;, d) 
/ \id 

From the definition of the w(A;, d) we see, that for k G (id, ic): w_|_(/c, d) = i|w+(A;, d)\. We make change of 
variables A: = iy in the last integrals and get 

where 

Ud)=j , ^^ 

y V(c2-y2)(rf2_y2)' 

and 

c 

Ay 



Io(d) 



^(C2_y2)(y2_rf2) 
d 



Yl 



Let make change of variable y = sd in the li{d). Then we see that li{d) is increasing function on d. 

Indeed, 

1 

Ii{d) 



6 J{l-s^)[l-Ks 



Now we make change of variable y = d+{c—d)s in the integral lo{d). Then we see that Io(d) is decreasing 
function on d. Indeed, 



Ud) 



ds 



y/{c{l + S) + d{l - S)){CS + d{2 - S))y/s{l - s) 



Then the statements of lemma is evident. 



n 



5.1 Proof of lemma 12.11 

Proof. To show that this system has a unique solution let us define the function 

1 



F(^,d) 



/ if^' - x'd^) 



1-A2 
2 - A^d^ 



dA. 



It is easy to see that there exists a function ^ = fi{d) such that F{fi{d),d) = and < fi{d) < d. 
Moreover, /x((i) is strictly increasing as d € [0, c]. Indeed, one can check that F{fi, d) is strictly increasing 
in /x and is strictly decreasing in d as < ;U < d < c. Now if < di < 1^2 < c then F {fj,(di) , di) = = 
F(/x((i2), ^2) < F {lJ'{d2) , di) , that is F{fi(di),di) < F{ii{d2),di) and, hence, /x(di) < ^{d2)- Furthermore 
jjL^d) is continuous function that follows from the representation: 



^?{d) 



x^d^ 



1-A2 

C2 - A2d2 



dA 



1-A2 



A2(i2 



dA, 



(5.27) 



which is equivalent to equality (|2.1U|) . Equation (|5.27p yields that ;u(0) = and ^(c) = —p. Hence 

v3 



d2 
^2(d) H is vary over segment 



0, 



5c^ 



as d is vary over segment [0, c] . So for any ^ G 



^ "1 
2 ' 3 



exists a single d G (0, c) such that (|2.10p and ()2.1ip are fulfilled. Let us denote this function by 



/: 



-c2 c2 



[0,c], 



/(O = d. 



-c2 c2 



there 



[0,c]. 



Moreover, / is increasing function. And therefore there exists an inverse function / ^ 
We can now rewrite equality ()2.1ip as follows: 

The last equality and the fact that ^{.) is increasing imply that /(.) is continuous function. So, system 
(j2. 1012. lip has a unique solution ^ = /x(/(C))i d = f{S,), and /.f(/(0))/(0 ^^^ continuous and strictly 
increasing functions. D 
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5.2 Asymptotic expansions of r*, Bg, A 

In this section we don't always write dependence of some functions on their argument to simphfy the text. 

Let us transform expression for A(d) (|2.14p : 

' nog(a+(A;)a„(fc))dfc pog{a+{k)a_{k))dk 

..n id w+{k,d) I w+{k,d) 12(d) ,__. 

^^^^ ^f^^T p^T" ~'- "^^ ^ ^ ^ 

w(A:,(i) w(iy,d) 

5.3 Expansion of I2 in r/ 

First we will treat l2{d). As for A; G (id, ic) w^{k,d) = i\w{k,d)\, then 

c 

W = |A: = iy|= /i^|MEtMii = |, = <i+(c-d).| = 

d 
1 

, ,, /" log(a+(i(d+(c-d)s))a_(i(d + (c-d)s)))ds 

(c-d) 







V^(c - (i)(l - s)(c + d + (c - d)s){c - d)s{2d + (c - (i)s) 



log (a+(i(d + (c - (i)s))a_(i(d + (c - d)s))) ds 



f log (Q- 

i A/i(l 



and 



v/s(1 - s) v^(c + d + (c - (i)s)(2d + (c - (i)s) 



1 
J /^/ ^.N ^ f log («+(i(c(l - ??) + C7?5))a_ (i(c(l - r?) + C7?g))) ds 

i Vs(l - s) v^(c(2 - ??) + C7?s)(2c(l - r/) + cr?s) 

1 
1 /■ log (a_|_(i(c(l — ?]) + C7ys))a_(i(c(l — fy) + c'i]s))) ds 



cj y/s{l - s)y/{2 -7] + r]s){2 -2r) + rjs) 



Now treat log a+a_. For k £ (ic, id) 



vri 



>c_{k) = |x(A;)|e4 



TTl 



x+(fc) = \>c{k)\e 4 
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(see the definition of >c(k) (??)) and then (see the definition of a{k) (??)) 

vri 

vri 



a.{k]a^{k]^U\x{kf + ^±-)=\k = iy\ = ^(.^+ '° " 



4 \ |x(A;)py 4 \\J c — y y c + y 



-w 1 H = — , =\y = d+ [c — d)s\ 

4.y c-y V c + yj 2^^^ ' V ; I 

|d = c(l-r?)| = 



1^{c - d)(l - s)(c + d + (c - d)s 



Then 



2^cr?(l - s)(c(2 - r/) + c??s) 2^/77(1 - s)(2 - r? + 77s) 

log (a_(^)a-.(^)) = log ,^^(,_^/(,_^^^^) = 

= "2 log??- 2^og(l - •s) - 2 ^°S(2- ^ + ^s) -log 2. 



Then the integral Y^idij])) is equal to 



I- —- logry — - log(l — s) — - log (2 — r/ + r/s) — log 2 



l2(d(r?)) = / -^ , ^ , ^ ds = 

J c^s{\ - s) v^(2 - 7/ + r/s)(2 - 2?? + ??s) 

— logr?ds I -— log(l-s)ds 



V^s(l - s) v^(2 - r? + 7?s)(2 - 2r? + r/s) 7 ^s(l - s) y/(2 - r? + r?s)(2 - 2?7 + r/s) 

/■ -77 log (2 -77 + 775) -log 2 

+ / r ds = 

J Ca/s(1 - s)^(2 - ?7 + r/s)(2 - 2r/ + r/s) 

1, ;G+o<"')-'^ 

= logr? / -^ ^ h 

i-^log(l-s)Q + 0(77))d. 

^ ^■log2 + 0(?7) 

-ds = 



Q cVs(l - s) 



-^logr7(| + 0(77))-l(-^log2 + 0(77)) + U-^^^ + 0(r7) 
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TT . TT log 2 „ , , TT , „ ^ / N 

= -^logr? ^^ + ^(^) = -^log27? + 0(r/). 

So, 

l2(d(7?)) = -f log27? + 0(r?). (5.29) 

4c 



The expansion of the Ii is more difficult than expansion of the I2 and is based on a ©-function iden- 
tity. 



5.4 Expansion oi T*{d{r])) in rj 

We know (by using Poisson summation formula), that 






47r2 

where r* = and from 1171 we know that 

r 



G(0|r((i)) c + d 



G(7ri|T(d)) Vc-d 

(see formula (4.34) in |17j). 

Let us recall that r](d) = 1 . Then, by using the Poisson summation formula (j3.15p . we get 

c 



2 -77(d) c + d e(0|r(d)) e{0\T*{d)) 



V r?(d) Vc-d e(vri|T(d)) Qfr^\r*(d))e^p^^' 

Now we use the inverse function h{.) of the function r* (see lemma ISTT]) and rewrite the last formula. 



l2-r]{h{T*)) e(0|r* 



Since 



"("(-•» e(^|.-)exp^ 



0(z|r*) = > exp < —T*m + zm > , 
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then 



viHr*)) = 2 



V I 



e^(0|r*) 



eMylrMexp^ 



+ 1 



E exp <j -T*m 



m=—oo 



*™2 



\-^ 



E exp <* -T*m{m +l)\) exp — 



+ 1 






1 + e^ 



*\\2 



+ 1 



8e"r +0 e~ 



4e-(l + 0(e^*)) 



Now we get 

T] = 8e^ + O (e^j , 

e^ = ^ + 0(r?2), 

r* = T*((i(??)) = -41og- + 0(7?). 



Also we have 



Vc2 - /i2(T*)e-i- = Vc2 - c2(l - r?(/i(r*)))2e-i- = 4c (l + O (e"^ 



5.5 Expansion of Ii in 77 

In the other way, (see the proof of lemma WT\i . 



T*{d) = -47r 



Ii(rf) 



and then 



Ii(d) = -^Io(rf). 



Now we get the asymptotic expansion of the Iq . 

c 

dy 



Ud) 



^(C2_y2)(y2_rf2) 



(c-d) 



47r 



I?/ = (i+ (c — d)s| 



ds 



(5.30) 



(5.31) 



(5.32) 



^(c -d){l- s){c + d+{c- d)s){c - d)s{2d + (c - d)s) 



ds 



^s(l - s) v^(c + d+{c- d)s){2d + (c - (i)s) 
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Let us make interchange of variables d = d{rj) = c(l — r?) in the last integral. Then 

1 

udm -- 



ds 



■\/s{l - s)y^(c(2 - 77) + C77s)(2c(l - r/) + crjs) 



ds 



cj ^s{l - s)y/{2 -T] + r]s){2 -2r] + r]s) c 



1 rQ + O^^))^" 
cJ ^/sil-s) 



2c 
And then by virtue of (|5.32|) and (|5.30p we conclude that 

hidiv)) = (^ + 0(^)) nog- + 0(7?)j =^log- + 0(r/logr/) 
Remark. Although in (|5.33|) we get only the first member of expansion of 



h{kc) = - [ ^ 



(5.33) 



, ^^^=^^^= as /c — )■ 1 (here k = —), but in this way we can get as much as desired 

V(l-A:2x2)(l-x2) ^ c^' 



members of this expansion. See also [28], problem 90. 



5.6 Expansion of A(d{rj)) 

Finally, by virtue of (|5.28p . (|5.33p and ()5.29p . we conclude, that 



and 



n 



A(d(r?)) 



-— log 2?7 + O (77) 
4c 



log - - log 2 + O (r?) 
-vr i-j 



I 



— log — h O (r/ log 7]) log — 1-3 log 2 + O (77 log 77) 

2c Tj 7] 



-vr 



V 



log- 



\ 



log2 77 



-Aidir^)) = -i 

TT 2 



i_ii^ + o^ 1 



log- 



log 77 



(5.34) 



In the following three paragraphes we get the expansion of /i. As we know from ()5.27p . 

Ad) = ° ^{^^^^ ^ 

a! , 7^2 _ ,,2 

/^^^=idy 
\/c2 - y2 
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And then 



J / o 9 ^ J / 9 9 ^ J / 9 9 -^ 



Ad) = ° y^^ - ° ^{^^ = d^ - ° ^(^^ =: d' - ^. (5.35) 

J I ry 9" y J / 9 9 ^ J / 9 9" ^ 







\/(? - y^ a/c^ - V^ \/c2 - y2 



5.7 Expansion of I3 in r; 

Let us first treat I3: 

d 



or 



T (A^ I V^IZa 
^3(^) = / Al 2 ^^ 



c{l-»?) 



I3(d(r?))= / ^^^S^^Zdy, (5.36) 

^ vc - y 

where d{i]) = c(l — r/). 

Let us differentiate /3((i(.)) in 77: 

(l3 0<(7?)= / . ^^^ ~'f ~^^ -^==dy = |see dOHDI = -c\l - v)h{d{v)) = 

J V *- (1 ~ 'n) ~ y \/c — y 

= |see (!5.33|) | = -c^(l - ??) ( — log - + (77 log?]) ) = — - log - + 0(r/logr/) . 

\2c ?7 J 2 7] 

Then 

hid{7])) =l3(d(0))+ / (I3 o (i)~(^)d^ = c - I / nog^ + 0(^log^)jd^=c-|r/log— + 0(r?2log??). 



(5.37) 

5.8 Expansion of L4 in r/ 

Now let us treat I4. 

Ud)-[^-^^^££^^y 

J v^^ 

^ ^ 

and 
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where d{r]) = c(l — r/). 

Let us differentiate l4((i(.)) in ij: 

c{l~r,) 



— — (r?) = / , „ dy= IseedOg) I = -3c (1-77)13 ((i(r/)) 

1^ J \ c^ — v^ 



„3n s /1 1 i„„8e , ^,2. 



I see ([071) 1 = -3c''(l-77) ( 1- -r?log h 0(77^ log?]) ) = -Ze' + O [rjlogr]) 



Then 



14(4??)) = 14(40))+ /'(Iod)~(^)d^=^c3-3c=^/' (1 + (^log^)d^= ^c^^-Sc^ry + O (772 logy?) . (5.38) 



5.9 Expansion of /i in 77 

Finally, by virtue of (|5.35|) and the fact that d = d{rj) = c(l — 7/), we get that 



.21, \ _ ,2 ^i{d) 



H\d) = d' 



Ud) 



and 



2 

-e' -3c^77 + (r/^ log 77) 
,\di,)) = c\l - ,f - -^^ 1 ^' ''' 

c - -77 log h O (77^ log 77) 

2 77 ' 

9 

2 , l--r/ + 0(77-log77) 

= C^(l-77)^ c^ ^ - 



So, 



and 



1 T] log h O (r/^ log 7/) 

2 7/ 

2/1 \2 2/-1 ^ 1 ^^ y^ / 2 1 2 \ 

c (1-^) - 3C {^-^r] + -r]log— + 0{r] log 77) 

c^(l - 77)^ + c^ ( -- + 37? - -7/log hO (7/^log^77) 

\ 3 3 77 

cM - + ?/- -??log h O (7?2 log2 r/) j =cM - - -77log— 2 +0(7/2log2?7) 

/i2(ci(77)) = c' Q - ^77log^ + O (T/^log^)") (5.39) 

3Ai2((i(?7)) / 8 ^/ 2i 2 nA 

= (l-77log-^ + 0(772log2r/)j. 



c2 



5.10 t; in 77 

Here we introduce one more small parameter v, which is more close in mean to ^ than 77. As we know 
from paragraph (|5.ip . 

r\d) = ^,\d) + ^-'^. 
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Let us define new variable 

v = v{^) = 1 - — , 

and remember the definition of the function d = /(^) from paragraph ()5.1 
Then 

and 

v{f-\dir^))) = 1 - ^ (^^'Kr/)) + ^ - I) = 1 - (l - r?log A + o (7?2log2,^)) - 

3 fd'^irj) c^\ , 8 3/ c?^(77)\ ^/ 2, 2 n 1 8 3 ,^ n ^/ 2, 9 x 
' ^'^ >-r/log^ + - 1 ^ +0 r^^log^r?) =7?log^ + -7?(2-r?) + r/2log2,? 



c2 V 2 27 '^7]e^ 2 V c2 7 ' ^' ^ '^ '^r]e^ 2 

1 8*2 r-v / 2 1 2 \ 

77 log h O (77 log 77) 

7/ ^ ^ 



or 



f = ^log- + 0(7?2log%). (5.40) 

8e 8e 7? ' 

5.11 ?7 in f 

Let us note that if x = x{y) is invertible function in some neighborhood of zero and 

x{y) = y + 0{y),y^O, 

then 



y{x) = X + O (x) , X — ;■ 0. 



V 77 8e 

We have x = — , y = — log — . Then 
8e 8e rj 



We have got a Lambert equation 

we^ = z, w <0,z <0, 

where 

w = log^, z = ^ + 0{v^). (5.41) 

8e 8e 

Following the [?], we get that this equation has the inverse for 77; << 0, and this inverse has the following 
expansion: 



L9 „ fL 



2 



where 



w = -Li-L2-j^ + 0{jj], z^-0, 



Li =log — , 

z 

L2 = log log — . 

z 
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Then 



qW =f^-LiQ-L2^^^ 



log 



|)exp(o(^M, . 

log log — / log^ log — 

1 ^ + O 



-0, 



log — 

z 



z 



z-f -0. 



And by virtue of (|5.4ip : 



8e 



log -Y, 



/ log log -Y, 



ii + o(.^) 



+ 



log 



V i, + 0(.^) 



log^ -Y, 



z^ -0, 



As 



v + Oiv^ 



( 



log- 



8e 



log log 



- + o(-^)^o(^°s'^°^^ 



+ 0(t;2) V ^v + 0{v'^) 



log- 



8e 



log' 



1 



t;^+0. 



then 



Be /8e 

^°g loS — , r-./ 2^ = 1°S log — (1 + O (u)) ) = log ( log — + O {v) 



log ( log — ( 1 + O . , 

V \ \ log V 



log log h O , , 

V \ log V 



log 



8e 



1 1 8e /, 2, 1 

log log — / log log ■ 

1 ^-^ + 0' ^ 



log- 



8e 



log' 



1 



, 1 I . logSe 
log - 1 + ^ 



log- 



, , , 1 / log8e 

log I log - I 1 + ^-^ 



log- 



V 



, 1 I , logSe 
log - 1 + 



\ 



log 



1 



V 



J 



I 1 o , , 1 log8e 
' log 8e + log log — I 



log- 

V 



( 



+ 



log log- 



1 



loe 



1 



1 



log- 

V 



^ log - / log^ log - 
" + O I f 



log'- 

V 



log- 

V 



1 / 2 1 

log 8e + log log - / log log - 
1 , ^ + ' ^ 



log- 

V 



log^- 

V 



v^+0, 
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and so we get that 



1 /. ., 1 



^ / log 8e + log log - / log^ log 

V = 



1 J ^ + 0| ^r I h ^^+0- (5-42) 



log- 

V \ V \ V 



\ log- \ log 

\ V \ 



5.12 Expansion of Bg 

As we know, (see (|2.12p ). 



Bg{d) = 24 I 



w{k, d) 
id 

We will get the asymptotic expansion of Bg as d is tend to c. 
Let us do interchange k = [y,y £ {d, c) in the last integral: 



Vc - y 



d 

1 



24(c -d) f 



{{d + (c - d)s)2 - /i2(rf))^(c _ d)s{2d + (c - (i)s)ds 



V^(c-d)(l-s)(c + d+(c-d)s) 



1 



2A(c-d)l ! ' iid + ic-d)s)^-f,\d)W{2d+{c-d)s) ^^ 







Now we make integrchange d = d[ri) = c(l — r?) and recall (see l5.39"|) . that 



o, „ .. c 



,2 



l^\d{ri)) = -{l + 0{ii\ogii)). 
Then 



B,,.(,), ^ 24C, / y^ V JU.-r,^,J 



1 

= 24c7?y XII 0c2 (1 + (r/logT?))") (1 + (r?)) ds = 



1 1 

= 24cr? I \-r^ f ^c^ (1 + O (j^log??)) j ds = IGc^? f J -^ {1 + O {i]logi])) ds = 



= 87rc^r/(l + 0(r?logr7)), r/ -> 0. 
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